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ABSTRACT 

We present an improved method for predicting the Sunyaev-Zeldovich (SZ) effect in 
galaxy clusters from spatially-resolved, spectroscopic X-ray data. Using the depro- 
jected electron density and t emperature profiles m easured within a fraction of the 
virial radius, and assuming a iNavarro. Frenk fc Wh ite (f995) mass model, we show 
how the pressure profile of the X-ray gas can be extrapolated to large radii, allow- 
ing the Comptonization parameter profile for the cluster to be predicted precisely. 
We apply our method to Chandra observations of three X-ray luminous, dynami- 
cally relaxed clusters with pubhshed SZ data: RX J1347. 5-1145, AbeU 1835 and AbeU 
478. Combining the predicted and observed SZ signals, we determine improved es- 
timates for the Hubble constant from each cluster and obtain a weighted mean of 
Hq = 69 ± 8kms~"'^ Mpc~^ for a cosmology with fJm = 0.3 and Q\ = 0.7. This result 
is in good agreement with independent findings from the Hubble Key Project and the 
combination of cosmic microwave background and galaxy cluster data. 

Key words: X-rays: galaxies: clusters ~ cosmic microwave background - cosmology: 
observations - distance scale - galaxies: clusters: individual (RX J1347. 5-1145, AbeU 
1835, Abell 478) 



1 INTRODUCTION 

' The inverse Compton scattering of cosmic microwave back- 
, ground (CMB) photons by hot electrons in galaxy clus- 
' ters leads to a distortion of the CMB spectrum along 
the line of sight, known as the Sunyaev-Zeldovich (SZ, 
ISunvaev fc Zeldovichlll972l) effect. The magnitude of the SZ 
effect is determined by the Comptonization parameter of the 
cluster gas, y{r), which is proportional to the line of sight 
integral of the gas pressure. 

It was recognized sw i ftly llSilk fc White! Il978t 
ICavaliere. Danese. fc de Zottil Il979t) that for an as- 
sumed, simple geometry, the combination of X-ray and SZ 
observations can be used to measure the angular diameter 
distance to a cluster. The ratio of the observed (based on 
radio/sub-mm observations) and predicted (based on X-ray 
observations) SZ signals is proportional to the square of 
the angular diameter distance, making this, potentially, an 
exceptionally powerful technique for probing the cosmic 
distance scale. 

Although the SZ effect is now empl oyed frequently 
to de t ermine extragalac t ic dis t ances (e.g.. iMauskopf et al, 
I2OOOI: ICarlstrom et al.l I2001I: iMason. Mvers fc Readheac 
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l200lt I Jones et alJ 1200 it iReese et al.l I2002H . most studies 
to date have relied on the app lication of the /3-model 
iCavaliere fc Fusco-Femiandll97a) in their X-ray analyses. 
The /3-model provides a simple analytical approximation to 
the spatial distribution of the cluster gas. These studies have 
also generally assumed isothermality and relied on broad 
beam measurements of the mean, emission-weighted tem- 
perature when calculating the gas pressure. 

With the Chandra X-ray Observatory and XMM- 
Newton it is now possible to measure the temperature, den- 
sity and pressure profiles of the cluster gas precisely. The sur- 
face brightness distribution, in particular, can be resolved on 
sub-arcsecond scales. Chandra observations have shown that 
/3-models do not, generally, provide good descriptions of the 
X-ray gas distributions in the cores of dynamically-relaxed 
clusters, and that the gas temperature is not isothermal but 
drops sharply within cluster cores (e.g ., ISchmidt et al .1200 ll : 
lAllen et allboOltlKaastra et alJl2004ll . The main limitation 
of the new X-ray observations is that (due to background 
levels and restricted fields of view; 8x8 arcmin'^ for the 
Chandra ACIS-S detector) the gas temperature can only be 
measured directly out to radii r < one third of the virial 
radius for most clusters. Since predictions of the SZ effect 
require the line-of-sight pressure integral through the clus- 
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ter, extrapolation of the X-ray results is therefore required 
for combined X-ray/SZ work. 

In this paper we present a recipe for calculating the pre- 
dicted SZ effect in clusters which makes use of the spatially- 
resolved spectroscopic techniques described in our earlier 
work, and which includes a new method for extrapolating 
the X-ray pressure profiles beyond the directly observed re- 
gion. Analytical formulae for the electron density and tem- 
perature profiles of the X-ray gas at large radii are obtained 
under the assumption of hydrostatic equilibrium. These can 
be computed easily and used for SZ analysis. We show that 
our method, which takes full account of the observed density 
and temperature profiles of the X-ray emitting gas, provides 
a significant improvement in accuracy with respect to the 
standard isothermal /3-model. 

The outline of this paper is as follows: In Section 2 
we review the standard /3-model approach and describe our 
improved method. In Section 3 we apply our method to 
Chandra X-ray and SZ observations of three highly X-ray 
luminous, dynamically relaxed galaxy clusters: RX J1347.5- 
1145 (z = 0.451), Abell 1835 {z = 0.2523) and Abell 478 
{z = 0.088). In Section 4 we use the combined data to 
measure the Hubble constant. In Section 5 we compare our 
results to those obtained using the standard isothermal /3- 
model approach. A summary of the results can be found in 
Sect. El 

Unless stated otherwise, all quantities are are quoted for 
a Hubble constant Ho = 70kms^^ Mpc~^, matter density 
Qm ~ 0.3 and vacuum energy density f^A ~ 0.7. Error bars 
correspond to 1 cr (68.3%) confidence. 



light, and rric the electron mass. Eq. ^ can b e integrated 
analytically for the isothermal /3-model (e.g., iBirkinshawl 
il999) to give 



y{R) = yo 



1 + 



with the normalization constant 
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jMauskopf et alJl2000h . where B{a,b) = r(a)r(6)/r(a + 6) 
is the beta function. 

In previous work, the temperature T has usually been 
taken to be the mean emission-weighted temperature of the 
X-ray gas, determined with broad-beam instruments. How- 
ever, as discussed in Section 1, observations with the Chan- 
dra and XMM-Newton satellites have shown that regular, 
dynamically relaxed galaxy clusters are not isothermal. This 
temperature variation should be accounted for in the anal- 
ysis.^ 

Finally, clusters are not infinitely large, as is usually as- 
sumed when calculating the predicted Comptonization pa- 
rameter profile using the /3-model. In what follows we show 
how modern X-ray data can be used to integrate eq. ^ 
without the assumption of isothermality, taking account of 
the finite size of galaxy clusters. In particular, where the re- 
gion of the cluster directly observed in X-rays is small, the 
use of an extrapolation procedure like the one described here 
can become important. 



2 PREDICTING THE SZ EFFECT 

2.1 The standard /3- model approach 

Within the context of the isothermal /3-model 
jCavaliere fc Fusco-Femianol Il978h . the surface bright- 
ness profile of a galaxy cluster can be written as 



I{e) oc 



1 + 



-3/3+i 



(1) 



where 9^ is the angular core radius and /3 is the slope pa- 
rameter. By calculating the X-ray e mission due to hot clus- 
ter gas at a tempe rature T (e.g., iKaastra fc Mewelll993t 
iLiedahl et al.ll^9^5^ ■ one can invert this to yield the cen- 
tral electron density rieo of the corresponding intrinsic (3- 
dimensional) electron density profile 



no(r) 



rico 



1 + 



-, -M 



(2) 



where the core radius Vc follows from the angular diameter 
distance d of the cluster, rc — d x 6c. 

The Comptonization parameter along a line of sight at 
an impact parameter R is defined by 



2 (Tt 



nc{r)T{r) 



r dr 



(3) 



(e.g., lBirkinshawl[l 999. where we have converted the equa- 
tion into an integral along radius r). Here ctt is the Thom- 
son cross-section, fee the Boltzmann constant, c the speed of 



2.2 An improved approach 

2.2.1 Discretization of the y-parameter calculation 

In order to calculate the expected y-parameter profile (eq.|3 
for a particular galaxy cluster, the temperature profile T{r) 
and the electron density profile nc{r) need to be known. In 
practise, d eprojection analyses of Chandra or XMM-Newton 
data (e.g.. ISchmidt et alJl200ll and references therein) pro- 
vide the electron density and X-ray gas temperature in 
discrete shells, with inner and outer radii rin.i and Vout.i, 
respectively. We can thus recast eq. Q into a sum starting 
with shell j at impact parameter R: 

yiR) 



nic c^ 

N 

+ ' 



out,j 
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E 
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For clusters with temperatures above about 8 keV, rel- 
ativistic corrections also become significant. Analytic 
and fitting relations for the relativistic corrections to 
the SZ effect have been worked ou t by several groups: 
|ltol^Kohjam^^^Nozaw^_^ |^9^; IChallinor fc Lasenbvl 
(199^; I Sazoiio v & Sunvaev (1998J. Using the analytical 
formulae of [ Challinor fc Lasenbvk (■199&'l up to second 
order, we determine the relativistic correction factor x 



Uncertainties in the determination of the distance scale as- 
sociated with the ass umption of isothermality were p revi- 
ously discussed b y, e.g ., Ilnagaki. Suginohara. &: Sutol ilQQST) and 
IVoshikawa ct al.i il998l) . 
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that needs to be applied to each summand of @, so that 



2.2.2 Extrapolation of the temperature and gas density 
profiles 

In the case of the Chandra X-ray Observatory, the detec- 
tor size and the particle background limit the regions of 
clusters for which direct temperature measurements can be 
made (r ^ ri) to a fraction of the virial radius. In order to 
calculate the predicted Comptonization parameter precisely, 
we need to extrapolate the temperature and electron density 
profiles out to the edge of the cluster, r2 - or at least past 
the point where significant contributions to y{r) are made. 
Here, for convenience, we will set ri = r25oo and 7-2 = r2oo 
(corresponding to the radii within which the mean enclosed 
mass density is A = 2500 and A — 200 times the critical 
density of the Universe pcrit(z) at redshift z, respectively). 
''2500 is typical of the outer radii for which useful information 
on the gas temperature profile can be obtained from Chan- 
dra observations. r2oo is used to mark the outer edge of the 
cluster. (In Sect.|S|we shall show that the precise choice of 
r2 does not affect the results significantly.) 

We parametrize the total mass distributions in the clus- 
ters using an NFW model 



p(r) = 



Pciit{z) Sc 



(7) 



(r/rs) (1 + r/rs) 

where p(r) is the mass density, r^ is the scale radius, pcrit = 
^^Iq , H{z) is the Hubble constant, G is the gravitational 



constant, Sc — ^ in{i+c)-c/(i+c) ' 



and c is the concentration 



parameter with c — r2oo /fs 

The temperature T{ri), electron density ndri), and gas 
mass-weighted temperature Tm.ri within ri are assumed to 
be known from direct spatially-resolved spectroscopy. We ex- 
trapolate the electron density from ri to r2 using a power- 
law model ne(r) ~ r~^. The hydrostatic equation is then 
used to determine the temperature solution T{r) for that 
power-law electron density profile in the given NFW poten- 
tial (see Appendix IXl. 

We fit the exponent 7 of the electron density profile so 
that the model complies with two fixed quantities: 

(i) the observed temperature T{ri) at ri 

(ii) the gas-mass weighted temperature Tin,r2 at ^2. 

The power-law model approximates the true electron den- 
sity profile between ri and r2. (It is assumed that the galaxy 
cluster ends at r2.) To determine th e gas-mass weighted tem - 
perature Tm.ra within r2, we use the lLokas fc MamonI (120011) 
solution for the radial dependence of the ratio of kinetic en- 
ergy Wkin and potential energy Wpot in an NFW potential 
with an isotropic velocity dispersion. The key assumption 
here is that the gas-mass weighted temperature Tm.r inside a 
radius r is proportional to the kin etic energy of th e mass dis- 
tribution inside this radius fe.g.. lEke et aljfl998l) . We then 
calculate the change, q, of the ratio VKkin/Wpot between ri 
and r-2 as 



q = 



W^kin(ri)/Wpot(r-i) T,„,.J [Mtot(ri)/ri] 



VKkin(r2)/Wpot(r2) T,^,rJ [Mtot (ra ) ^2] ' 



(8) 



Table 1. Details of the Chandra NFW mass models 
Object redshift c Tb (Mpc) 



RX J1347.5-1145 
Abell 1835 
Abell478 



0.451 

0.2523 

0.088 



6.34; 
4.21 

3.88 



1.61 
1.36 

+0.53 
0.61 
0.28 
-0.36 



^•■^'-0.12 

'i^+o-i* 

^■^^-0.09 
1+0.12 
0.07 



0.61'' 



quantity q can be computed as a function of the scale radius 
rs and the concentration parameter c of the NFW profile 
using the analytical formulae included in Appendix El 

Thus, given both the gas-mass weighted temperature 
rm.ri at ri and a specific NFW mass model from the Chan- 
dra data, the gas-mass weighted temperature Tm.r at r2 can 
be obtained from eq. (|5J. For the correct electron density 
slope 7, this is equal to the temperature that follows from 
the temperature solution (eq. IA3t and gas density profile 
Pg(r) (for which we assume Pg(r) — 1.1345 rrip nc(r)) via 
the integral 



Mg,,,T,„,,i +A7V / T{r)pg{rydr , (9) 



where Mtot{r) is the total enclosed mass within radius r. The 



where Mg^r is the gas mass inside the radius r. 

We note that in cases where the X-ray data only extend 
to radii ro <C »'25oo (within which cooling and/or heating 
effects may have modified Tm.ro significantly), one should 
extrapolate the gas density (Appendix ^ from ro to r2oo, 
use this extrapolation to estimate the gas temperature and 
density at ri = 7-2500, and then apply eq. © between r25oo 
and r2oo as usual. (This is possible because the extrapola- 
tion recipe is attached continuously to the Chandra data 
so that Tm.r can be calculated for any exponent 7 of the 
electron density extrapolation, regardless of where we at- 
tach the extrapolation.) We show below that this approach 
leads to robust answers in the case of Chandra observations 
of Abell 478, where the data cover only a relatively small 
radial range. 



3 APPLICATION TO PUBLISHED CHANDRA 
AND SZ OBSERVATIONS 

In this section we apply our method to Chandra and 
SZ observations of three galaxy clusters: RXJ1347. 5-1145, 
Abell 1835 a nd Abell 47 8 . The X-ray data were origi nally 
pubhshed bvlAllen et al.1 J2002h . ISchimidt et al.1 J200lh and 
ISun et al] i2003l) . respectivelv. 



3.1 Analysis of the X-ray data 

Inside the region accessible to Chandra, the deprojected 
temperature and electron density profiles for the X-ray gas 
can be determined directly (under the assumption of spher- 
ical sy mme try) using the metho ds described by Allen et sil\ 
hOOW and lSchmidt et al.l J200lll . These data are then used 
to calculate the 68 per cent (la) confidence region of NFW 
mass models that provide the best-fit to the Chandra data. 
The resuhs on the NFW mass models for RX J1347.5-1145, 
Abell 1835 and Abell 478 are summarized in Tab.0 

We have carried out the extrapolation procedure de- 
scribed in Sect. 12.231 for every mass model within the 68% 
confidence region obtained for each cluster. The radii ri and 
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Figure 1. Extrapolated temperature and electron density profiles for RX J1347. 5-1145, Abell 1835 and Abell 478. The best fit profiles 
are indicated by the solid lines. The upper and lower envelopes of the model ranges are plotted with dashed lines. For clarity we have 
only plotted the profiles up to the smallest A = 200 overdensity radius of each model ensemble (see Tab. 1^. 



Table 2. Details of the cluster model extrapolation 

Object ri r2 T(ri) Tm.ra nc(ri) 

(Mpc) (Mpc) (keV) (keV) 10" 



-3 



RX J1347.5-1145 O.JStrol ^■^'^tlH IS-lSlJ-^? dSdtlH ^AAt^ofo l-25t^« S.SOt? ^^ 
Abell 1835 0.67tlf, 2.36«;?^ 11.34«;S 7.97+_r,l 7.461^4^ 1.361^^ 1.88l°;ro 

Abell 478 2.36«:1? 8.29l°:« 6.41«;« 5.92t»:» l-^Ol^^ l-90lH^ 



r2, the temperat ure T^, , the ratio q (e g. El calculated using 
the formulae in iLokas MamonI i2001il and the effective 
power-law exponent 7 of the electron density distribution 
for the clusters are given in Tab. H In the case of Abell 
478, the Chandra data do not extend beyond r ~ 0.3h^Q 
Mpc (corresponding to an overdensity A = 8500). Within 
this radius radiative cooling may have affected the observed 
temperature of the cluster gas. Following the instructions 
from Sect. l2.2"^ for such a case, the temperature profile and 
electron density profiles were extrapolated from rgsoo, but 
the exponent 7 was determined using eq. ||HJ between r25oo 
and r2oo. 

The application of our extrapolation procedure leads 
to the temperature and electron density profiles shown in 
Fig. In this figure, both the best-fit profile and allowed 
range of models are shown. Note that the small scatter in 



the electron density profiles is a consequence of the small 
scatter in the luminosity profile derived from the Chandra 
data. The temperature profiles have been binned using a 
simple emission- weighting scheme. 

For Abell 478, it is also possible to determine the power- 
law slope of the electr on density di stribution beyond ri using 
ROSAT observations llAllenl2000l) . Fitting the ROSAT data 
between 0.3 Mpc and 0.9 Mpc yields a slope of 7 = 1.84 ± 
0.07. This result is consistent with our model predictions 
given the expected steepening of the profile beyond the limit 
of the ROSAT data. 

Finally, we note that the power-law model used in ex- 
trapolating the density profile provides only an approxima- 
tion to the density values at large radii. (The true profile will 
steepen with increasing radius.) However, the pressure pro- 
file - the relevant quantity for predicting the Comptoniza- 
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tion parameter profile - should be predicted accurately by 
the extrapolation procedure over the full range of radii. 

3.2 Combination with the SZ data 

3.2.1 RX JI347.5-II45 

Detections of the SZ effect in R X J1 347.5-1145 were 
pubhs hed bv iPointecouteau et all ll200ll POl) at 142.9 
GHz iKomatsu et all toOlTl at 21 GHz and 150 GHz 
and iReese et al.l l2002'. R02) at 30GHz using the Diabolo, 
Nobeyama and Owens Valley Radio Observatory (OVRO) 
instruments, respectively. The Comptonization parameters 
of all groups are consistent with each other. We compare the 
Chandra prediction with the Diabolo detection, which has 
the smaller beam size at 142.9 GHz of the two bolometers 
(Diabolo and Nobeyama), as well as a smaller error bar on 
the central Comptonization parameter than the OVRO in- 
terferometer result. POl also published the flux decrement 
of the background radiation in four quadrants, which al- 
lows us to exclude the south-east quadrant, which contains 
hot, probably shock ed, gas (Komats u et al . 2001; Allen et 
al. 2002). Recently, iKitavama et alJ (|20o3|) have also car- 
ried out a detailed study of the intracluster medium in RX 
J1347.5-1145 using SZ observations. They obtain results in 
good agreement with those determined from the Chandra 
X-ray data. 

In order to avoid the complicating effects of the shocked 
gas to the south-east of the cluster centre in RXJ1347-1145, 
only data from the other three quadrants were used to de- 
termine the X-ray mass model llAllen et alJl2002ll . POl pub- 
lished the flux decrements in the four quadrants, exclud- 
ing a small 16.2 arcsec x 16.2 arcsec square region around 
the nucleus (to avoid contamination by the central radio 
source). Adding all four quadrants, they find fiux decre- 
ment of —(39.0 ± 2.8) mjy in a square 100 arcmin x 100 
arcmin area around the cluster centre. However, this be- 
comes — (26.8 ± 2.4) mJy once the south-east quadrant is 
excluded. This is the measurement we use for comparison 
with the X-ray data. 

The intensity change Al due to the SZ effect can be 
calculated from the Chandra model. In the Kompaneets ap- 
proximation (e.g., .Birkinshaw..l999.'l this is given by 



Al = loy g{x) 
where Jo = 2(fcB Tof/h^c^, 
.„ . "a;(e" + l) 

9{x) = 



(10) 



(11) 



and X — Iw /k^To, with the Planck constant h, observing fre- 
quency f and microwave background temperature To . Where 
relativistic corrections are significant, the Comptonization 
parameter y has to be modified appropriately (Sect. l?T2.H . 

3.2.2 Abell 1835 

T he detection of the SZ e ffect in Abell 1835 was published 
bv lMauskopf et al.l J200l MOO) using the Suzie I and Suzie 
II bolometers. R02 have also published a detection using the 
OVRO interferometer at 30 GHz, which provides a consis- 
tent, but more precise determination of the Comptonization 
parameter. Here we compare our Chandra comptonization 



profiles with the co-added SuZIE I scans published by MOO, 
using their (j> = —0.1 arcsec offset of the X-ray cluster centre 
from the scan centre. To calculate their difference channels 
D3 (two beams separated by 4.6 arcmin) and T123 (triple 
beam chop of three beams separated by 2.3 arcmin) of the 
bolometer array we followed the description given by MOO. 



3.2.3 Abell 478 

AbeU 478 was obse rved bv iMvers et al.| Jl997ft at 32 GHz 
and reanalysed by iMason et al.l J200ll MMR). The ob- 
served decrement was published as an average y-parameter 
y = 7.52 ± 0.56 x 10""" within the telescope beam. Note 
that we use their value without the relativistic correction, 
as we apply this correction to the X-ray model. We compare 
the observed SZ decrement with the Chandra p rediction us- 
ing th e beam switching technique described bv lMvers et alJ 



4 HUBBLE CONSTANT DETERMINATION 

The Comptonization parameter y defined in eq. Q depends 
upon the sq uare root of the angular diameter distan ce to 
the cluster dSilk fc White! ll978l:ICavaliere et al]|l979H . By 
comparing the observed Comptonization parameter, j/obs, 
within a given beam or aperture with the predicted value, 
J/prod, from the X-ray data (for a given cosmology with 
Ho = 70kms-^ Mpc"^), we can measure the Hubble con- 
stant e.g. (MMR), 



Ho = 



yprcd 
Vohs 



X 70kms"^ Mpc"^ 



(12) 



Since eq. 1)10^ is linear in the Comptonization parameter, 
these considerations also apply to the flux decrement AF^ = 
ujAIv fSect. IT^ . where u) is the solid angle of the emitting 
area, and AI^ is the intensity decrement. 

In order to compare our X-ray-predicted Comptoniza- 
tion parameter profiles with the direct SZ observations, we 
have convolved the X-ray profiles with the instrument beams 
for RX J1347.5-1145 (22 arcsec FWHM beam) and Abell 478 
(7.35 arcmin FWH M beam and 22.16 arcmin beam switch- 
ing as described bv lMvers et alJ lll997l) 'l. For Abell 1835 we 
have added an additional normalisation parameter to the 
NFW parameter space of scale radius and concentration pa- 
rameter constrained by the Chandra data (Table 0. To de- 
termine the best fit value and la error for this normalisation 
parameter we added the contribution from the SuZIE I 
co-added scans (1.7 arcmin FWHM beam and 2.3 arcmin or 
4.6 arcmin beam separation) to the contribution from the 
Chandra data. The normalisation parameter directly mea- 
sures the distance scale and is represented here by the re- 
sulting Hubble constant. 

The results for the three galaxy clusters are shown 
in Table |3 For quantities with two-sided error bars, the 
root-mean-square error was used to determine the error 
on the Hubble constant. The weighted mean of the three 
Hubble constant determinations in Table |3 is Ho = 69 ± 
8kms"^ Mpc"^ 
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Table 3. Results: Column two lists the observed quantity used to measure the Hubble constant. Columns three and four contain the 
measured and predicted values. The fifth column lists the observed areas. Column six contains the implied Hubble constant {Qm = 0.3, 
Ha = 0.7). 

Object quantity observed predicted area Ho [kms~^ Mpc~^] 

RX J1347.5-1145 AF^[mJy] -26.8±2.4 mjy -25.5 ± 2.3 mjy 3 quadrants 63.4 ±16.1 

Abell 1835 SuZIE I scan direct fit 24 arcmin scan 77.5 ±16.5 

Abell 478 5 [10-^] 7.52 ± 0.56 7.43+^'j4 7.35' FWHM beam 68.3 ± 10.5 



Table 4. /3-model fits to the Chandra data 

Object T Be rico l3 

(keV) (arcsec) (cm^^) 

RX J1347.5-1145 12.2 4.8 0.207 0.57 

Abell 1835 8.1 9.0 0.100 0.61 

Abell 478 6.8 26.6 0.068 0.55 



5 COMPARISON WITH THE /3-MODEL AND 
PREVIOUS STUDIES 

As the method for calculating the predicted Comptoniza- 
tion profiles presented here is significantly different from the 
standard isothermal /3-model approach, it is instructive to 
compare the results on y{r) and the Hubble constant ob- 
tained with the two approaches. 

We have fitted the Chandra surface brightness profiles 
with a /3-model (eq. inside a radius of 2.5 arcmin for RX 
J1347.5-1145, 3 arcmin for Abell 1835, and between radii of 
1.5 arcmin and 4.5 arcmin for Abell 478. Note that in the 
case of Abell 478 the core of the cluster had to be excluded, 
because the central surface brightness profile is not fiat, as 
is required by the /3-model. The mean, emission-weighted 
temperatures from single-temperature fits to the Chandra 
spectra and the /3-model parameters are given in Tab. |1| 
The temp eratures were determined using the MEKAL 
iKaastra fc Mewe 1993; Liedahl ct al. 199^) plasma model. 
The gas density normalization rieo was calculated from the 
surface brightness using the MEKAL plasma emission model 
with the spectroscopically determined metallicities. Rela- 
tivistic corrections were applied as described in Sect. 12.2.11 
Using eq. ^ , together with the appropriate relativistic cor- 
rections factors, we obtain the central Comptonization pa- 
rameters yo = 1.13 X 10"^ for RX J1347.5-1145, yo = 
4.2 X 10"^ for Abell 1835 and yo = 3.5 x 10"'' for AbeU 
478. 

Fig. |5| shows the best-fit Comptonization parameter 
profiles determined with our new method (solid line) to- 
gether with the results from the standard isothermal /3- 
model approach (dashed line). This plot shows that the 
Comptonization parameter profiles for the /3-model differ 
substantially from the profiles obtained using the new, de- 
projected/extrapolated solution. Firstly, the overall normal- 
ization in the core can be quite different. Secondly, the 
slope of the /3-model is too shallow at large radii due to 
the infinite extent assumed in that model. (We note that 
the /3-model can yield significantly different answers when 
applied to other data sets c overing different ranges of radii 
frointec outeau et al.lf200ll: iReese et all 120021: iMason et all 
L2OOI). Here we have applied the two approaches to the same 
Chandra data simply to enable a direct comparison of the 
methods on the basis of the same, well-defined data sets.) 

Applying the procedure described in Sect. 21 to the 




10 100 1000 

radius (kpc) 



Figure 2. Best-fit fully deprojected and extrapolated Chandra 
Comptonization parameter profiles (solid lines) . The results from 
/3-model fits to the Chandra data are plotted with dashed lines. 
All profiles are calculated for Ho = 70kms~^ Mpc~^, Cm = 0.3 
and Qa = 0.7. 



measured /3-profile parameters yields Hubble constant es- 
timates of Ho = 32kms-^ Mpc"^ for RX J1347.5- 
1145, Ho = 29kms-^ Mpc"^ for Abell 1835 and Ho = 
103 kms"^ Mpc"^ for Abell 478. One can immediately see 
that the application of the isothermal /3-model does not yield 
a consistent estimate when applied naively to Chandra data 
for the present sample of clusters. 

It is also interesting to compare these estimates with 
previously published /3-model estimates for the clusters dis- 
cussed here. We have translated the electron densities from 
other studies into our assumed cosmology (Q^n ~ 0.3, Qa ~ 
0.7 and Ho = 70kms-i Mpc'^). 

RX J1347.5-1145: POl determined a Hubble constant 
Ho = 44±6kms"^ Mpc"^ using a kT = 9.3 keV isothermal 
model with the /3-model parameters 9c = 8.4 arcsec, = 
0.56 and n^o = 0.102 cm"^ from ROSAT llSchindler et alJ 
^93). Using the same temperature and a slightly different 
/3-model, R02 found an angular diameter distance that im- 
plies Ho = kms~^ Mpc~^. This estimate is consis- 
tent with the value we obtain using our fully deprojected 
and extrapolated cluster model (Tab.j^J. However, the gas 
temperature and the other /3-model parameters used by R02 
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are not consistent with the /?-model parameters from Chan- 
dra (Tab.BJ. 

Abell 1835: For this cluster, MOO measured a Hubble 
constant Hq = 66l^^kms"^ Mpc"^ for a kT = 9.8 keV 
isothermal model with the /3-model parameters 9c — 13.2 
arcsec, /3 = 0.58 and noo = 0.058 cm~^ from a fit to ROSAT 
data. This temperature was obtained by allowing for a cool- 
ing flow component in the core of the cluster. R02 obtained 
an angular diameter distance to the cluster which implies 
Ho = 55.41^^8^ km s"^ Mpc'^ , based on a kT = 8.21 keV 
isothermal model (no cooling flow correction). The differ- 
ence between these two measurements can be attributed to 
the different isothermal temperatures and small differences 
in the assumed /3-model parameters. Note also that the /3- 
model parameters used in these studies are significantly dif- 
ferent from those determined from Chandra (Tab. ^ over 
a smaller range of radii. This explains the difference in the 
predicted Hubble constant. However, the MOO result is con- 
sistent with the result we obtain from the deprojected and 
extrapolated X-ray data, which shows that it is possible to 
partially correct the /3-model in the central region of the 
cluster when the surface brightness profile observed with 
the ROSAT field of view is used. For larger radii this ap- 
proximation will fail, however, because of the fixed /3-model 
slope. 

Abell 478: MMR measured a Hubble constant Ho = 
6II21 kms~^ Mpc~^ assuming an isothermal kT — 8.4 keV 
model with the /3-model parameters 9c = 1 arcmin, /3 = 0.64 
and Uco ~ 0.023 cm~^ from ROSAT. This temperature also 
acco unts for th e pres ence of a cooling flow in the cluster 
core. ISun et aTl i2003h improved upon the result by MMR 
by determining the Comptonization parameter as a function 
of radius from a deprojection of the Chandra data, and by 
extrapolating the Chandra data with a /3-model based on a 
fit to the combined surface brightness profile from Chandra 
and ROSAT. They used /3 = 0.68 and obtained a Hubble 
constant estimate Ho = 64l5g kms~^ Mpc~^, which lies be- 
tween our value and the one obtained by MMR. 



6 SUMMARY AND DISCUSSION 

We have presented a new method for extrapolating the pres- 
sure profiles of galaxy clusters beyond the region accessi- 
ble to current X-ray satellites such as Chandra and XMM- 
Newton, for use in combined X-ray and SZ studies. Our 
method assumes hydrostatic equilibrium, an NFW mass 
profile for the clusters and predicts the temperature pro- 
file and an effective power-law slope of the electron den- 
sity profile. The method does not assume isothermality and 
provides a simple extrapolation recipe that can be easily 
implemented. We have applied our method to the Chandra 
data for three X-ray bright, dynamically relaxed galaxy clus- 
ters, RX J1347.5-1145, AbeU 1835 and AbeU 478, and have 
obtained detailed temperature, electron density and Comp- 
tonization parameter profiles for these clusters out to the 
virial radius. 

Since all three clusters have published SZ detections, 
we were also able to measure the Hubble constant for each, 
as shown in Table |21 The individual measurements are con- 
sistent with each other and yield a weighted mean of Ho = 
69±8kms~^ Mpc~^. This result is consistent with findings 



from the Hubble Key Project [Ho = 72±8kms"^ Mpc"^ 
iFreedman et alj I200J) and combined studies of the cos- 
mic microwave background and X-ray galaxy cluster data 
(Hq = 68 Atf° kms-i Mpc'^ lAllen. Schmidt fc Bridlj 
l2002). Note that our weighted mean value does not include 
additional systematical errors, due to, e.g., clumping or as- 
phericity, although for such regular, relaxed clusters these 
should be at a minimum. 

The contribution to the predicted Comptonization pa- 
rameter from material beyond the regions directly observed 
with Chandra, modelled by our extrapolation procedure, 
varies from system to system. For RX J1347. 5-1145 the cor- 
rection to the Comptonization parameter in the region used 
in Table El amounts only to 1 per cent and is negligible. For 
Abell 1835, however, the contribution from material in the 
extrapolated region amounts to 20 per cent, and for Abell 
478 the contribution rises to 68 per cent. This shows that a 
robust extrapolation method is necessary when combining 
X-ray and SZ observations to infer cosmological information. 

We have estimated the effect of changing the outer ra- 
dius r2oo by repeating the calculation for Abell 478, but 
stopping at r^oo- This leads to a 4.2 per cent smaller av- 
erage Comptonization parameter, well within the observa- 
tional error bar. 

We have compared the Comptonization parameter pro- 
files predicted by our method with the profiles obtained by 
fitting an isothermal /3-model to the Chandra data. We find 
that a naive application of the /3-model approach leads to 
profiles with significantly different central Comptonization 
parameters and predicts different slopes for the Comptoniza- 
tion parameter profiles at large radii. The isothermal /3- 
model also leads to inconsistent Hubble constant estimates 
when applied to the Chandra data sets considered here. 

It is clear that for galaxy clusters like Abell 478, a field 
of view larger than the one afforded by the Chandra ACIS- 
S3 detector would be beneficial. In this respect the com- 
bination of Chandra and XMM-Newton observations will 
help {e.g. Pointecouteau et al. 2004), although the detector 
background will still prohibit precise measurements of the 
temperature profile beyond r ~ 0.5r2oo. (We also note that 
the high spatial resolution of Chandra is important in re- 
solving the central temperature profile and constraining the 
best-fit NFW mass models.) 

Observations for a larger sample of clusters, as well as 
deeper X-ray and SZ observations, should make this a power- 
ful method to explore the cosmological distance scale. In the 
first case, it will be important to concentrate such studies on 
the largest, dynamically relaxed clusters for which system- 
atic uncertainties associated with the deprojection method 
and assumption of hydrostatic equilibrium are at a mini- 
mum. 
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APPENDIX A: AN EXTRAPOLATION 
SOLUTION FOR THE X-RAY GAS 
TEMPERATURE PROFILE IN AN NEW 
POTENTIAL 

Let the mass distribution of a galaxy cluster be parametrized 
by a NEW mass model (eq.Ql. Also, let the electron density 
distribution for a certain range of radii be parametrized by 
a power-law model 



(Al) 



The temperature as a function of radius can then be deter- 
mined under the assumption of hydrostatic equilibrium by 
solving the hydrostatic equation 



1 d(no kT) 



GM 



(A2) 



where k is the Boltzmann constant, jj, is the molecular weight 
of the gas and nip is the proton mass (note that the normal- 
ization of lAlll drops out of this equation) . Using the method 
of the variation of the constant one finds that is possible to 
write the solution of IIA2t in the form 



fcT(r) 



^ri/ \ ^ (l + 7)riri+T 



2Fi(-7,l;l-7; ) 



ri ^ 



-rir"+^2Fi(-7,l;l -7; -) 

-r,r\+-' ln(l + ^) + r,r^+-< \n(l + ^)))v3) 

with the boundary con dition T(ri) = Tj. 2F1 is the hyperge- 
ometric function (e.g., iGradshtevn fc Rvzhikl200(l . p. 995) 
which is well-defined on the whole negative real axis and is 
available in standard mathematical packages. 



APPENDIX B: THE ENERGY DISTRIBUTION 
OF THE NEW MASS DISTRIBUTION 



In this Appendix we list the lLokas & MamonI ||200J) expres- 
sions for the potential energy Wpot and the kinetic energy 
M^kin in an NEW potential with an isotropic velocity disper- 
sion, as a function of radius r. The ratio of these quantities is 
used in Sect. l2.2"^ to determine the cluster model extrapola- 
tion. It can be calculated from the concentration parameter 
c and the scale radius r^. 

The radial dependence of the potential energy asso- 
ciated with an NEW mass distribution with virial radius 
7'200 = c X Ts as a function of (scaled) radius s = is 



W^pot(s) 



2 ln(l + cs) 



1 + cs 



(Bl) 



(1 +CS)2 

(eq. 21 in lLokas fc MamonllioOlh . where Woo is the asymp- 
totic value 



Wo. 



GMl,{r 200) 



2r41n(l-Fc)-c/(l-Hc)]2' 



(B2) 
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The radial dependence of the kinetic energy, assuming 
isotropic velocity dispersion, is 

M/kin(s) = ivyoc{-3 + — ^ 21n(l + cs) 

A 1 + cs 

+cs[5 + 3 ln(l + cs)] - c^s^[7 + 6 ln(l + cs)] 

+c^s'^[n^ — In c — Ins + ln(l + cs) 

+3 In^ (1 + cs) +6Li2 (-cs)]} (B3) 

(eq. 24 in lLokas fc MamonlEoOlh . where Li2 is the diloga- 
rithm. 

This paper has been typeset from a Tff)i/ Wlft^K. file prepared 
by the author. 
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